Average size of random polygons with fixed knot topology 
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Abstract 

We have evaluated by numerical simulation the average size Rk of random polygons of 
fixed knot topology K — 0,3i,3iJ4i, and we have confirmed the scaling law R 2 K ~ N 2uk 
for the number TV of polygonal nodes in a wide range; N — 100 - 2200. The best fit gives 
2vk — 1.11 - 1.16 with good fitting curves in the whole range of N. The estimate of 2vk is 
consistent with the exponent of self-avoiding polygons. In a limited range of N (N > 600), 
however, we have another fit with 2vk — 1-01 - 1.07, which is close to the exponent of 
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random polygons. 



1 Introduction 

Topology of polymers is an important issue in understanding the physical properties of polymer 
materials, such as the viscoelasticity of polymer solutions [Q. It is also relevant to fundamental 
problems in biology |2j; the enzymes can cut and reconnect DNA strands, which are crucial 
functions in the mechanism of the complex processes of reproduction, transcription and recom- 
bination of DNA strands 1 3 . The topology of DNA rings can be used as a probe to detect the 
action of the enzymes on the DNA strands by observing the change of their knot type. The re- 
cent experimental observation [I] reports that the average polymer size depends on the topology 
of polymers. 

One important quantity in the study of ring polymers is their average size with fixed knot 
type. This quantity has been estimated in several numerical methods |S1 1^1 FTl IHIIUI fTUl fTTl fT^j . 
The renormalization group (RG) argument leads to the power law scaling for the average size 
of linear polymers as thier length increases |13[I14[ IT5] , 



where N is the number of segments, and R 2 (N) is the mean square size of polymers. The same 
scaling law should hold for ring polymers without topological constraint with the same value of 
the exponent v. It is commonly accepted that the size of ring polymers with fixed knot topology 
K is given by the same scaling relation as (fl|). 



This conjecture is supported by numerical simulations 7, 8, 9, TUlllllll2| . 



The scaling law of the form Q has already been discussed in connection with random 
walks (RW) and self-avoiding walks (SAW). Closed paths of RW and SAW are respectively 
random polygons and self-avoiding polygons (SAP). The exponents in these two cases are well 
established: vrw = \ and vsaw — 0.588 [TH1[T7|. The central question regarding the scaling (J2J) 




(1) 
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is: What is the value of uk ? Is vk related to either vrw or vsaw"? Note that the exponent 
i>K, or the coefficient Ar or Bk, should depend on the knot type K. 

It has been found by numerical analysis that, in the case of SAP, the exponent vk does not 
depend on the knot type K, and is given by vsaw [HIE]- The question is then whether the 
effect of topology appears in Ak or at the level of correction to the scaling (j2J, i.e. in Bk |18j . 
Simulations by van Rensburg et al. j7[|S] support the conjecture that knot topology does not 
affect the coefficient Ar. 

For the Gaussian model of polymers, des Cloizeaux made the conjecture that the constraint 
of topology should lead effectively to the growth of their average size, so called the "topological 
excluded volume effect" , even though it has no excluded volume interactions ^5] • It is predicted 
that the exponent for polymers with trivial knot topology obeys v§ < i>saw- Deutsch has 
performed a simulation for phantom chains in the case of trivial knot and obtained the value 
of f0 consistent with this inequality 0. The conjecture is also supported by a phenomenological 
argument j2H] for and other knots. 

Recently, Shimamura and Deguchi have evaluated the exponent vk for the radius of gyration 
for knot types 0, 3i and 4i in the case of the Gaussian random polygons 11 . Taking into account 
the correction term in the scaling P|). they have found that the fit gives vk — i'rw- 

In this paper, we re-examine the radius of gyration of the closed phantom chain model. 
For the topological effect on the size of random polygons, the two different answers have been 
presented. Furthermore, other possibilities are not excluded. The topological effect should also 
be investigated in different models. We employ the pivot algorithm which is a modification of 
the original one |5| so that the possible bias is reduced. The simulation procedure used in [S], 
which is a kind of pivot algorithm |211I22| for continuum models, tends to pick up extended 
conformation more frequently. It is conceivable that this bias leads to an artificial expansion of 
chains and to raise the exponent vk from vrw- We have evaluated the radii of gyration for the 
range of chains N = 100 - 2200 for three knot types, trivial knot, trefoil knot and composite 
knot, denoted by 0, 3i and 3i|4i. In the case of 0, our improved algorithm has given a result 
consistent with that of [S]. We have further results in other cases. 

We have found two good fits to our simulation data. The best fit gives 2vk — 1-H ~ 1-16. 
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It is consistent with vsaw^ which also agrees with the result of Deutsch: v§ ~ 1.17 [SJ. The 
second fit gives 2vk — 1.01 - 1.07, which is in accordance with vrw- The first fit gives good 
fitting curves for the whole range of TV" investigated in simulation: N = 100 - 2200. 



We consider closed phantom chains in three dimensions, consisting of N line segments of length 
a. We call them iV-node polygons. We assume that the polygons have no excluded volume. 
It means that the nodes and segments of the polygon are purely geometrical points and line 
segments, respectively (Hj. This model may be regarded as a model for polymers at the #-point 
or a polymer in a melt. 

A polygon Vn is defined by the set of position vectors of its nodes, Vn = (Ri, R2, ■ ■ ■ > Rn)- 
All cyclic permutations of one position vectors correspond to the same polygon. The vectors 
satisfy the geometrical constraint \Ri+\ — Ri\ = a for 1 < i < N (Rn+i = Ri)- Each polygon 
is topologically equivalent to a knot K in three dimensions. The configuration space C of the 
model is divided into subspaces Ck in which all polygons have a fixed knot type K, C = ^2 K Ck- 

The radius of gyration of polygons is given by 



This definition is different from that of the polygon size employed in |Hj but their asymptotic 
behaviors in the limit of large iV are the same. We set the segment length a to unity; macroscopic 
properties are independent of the microscopic parameter a. 

We generate a large number M of polygons with length N. They are used to evaluate the 
mean square size of polygons without topological constraint 



2 Model 




(3) 




(4) 



and the same quantities for fixed knot topology, 
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Here the indicator function x(^) K) takes the value 1 if V G Ck and zero otherwise. The number 
of polygons Mk with fixed topology K is Mk = Yli X^Pn,ii K). 
We are concerned with the scaling relation 

= ^f N ^K U + ABkN -i/2 + OiN- 1 )} , (6) 

where Avk = vk — v, ABk = Bk — B. It follows from the scaling laws (^Q) and © with 
A = Ak = 1/2 (BJ. We perform the fit simulation data using © 



3 Simulation Procedure 

Our sampling of the polygons follows the dynamic Monte Carlo method using the pivot algorithm 
|2HI22| . applied to the continuum model by Deutsch In the continuum model, a pivot move 
for a polygon is a rotation of a chain of segments randomly chosen from the polygon around the 
axis passing the two endmost nodes of the chain by a random amount of angle 9. 

In practice the method used in [§] imposes two restrictions on pivot moves, 1) the banning 
of self-intersections during the move and 2) the range of angle 9. There is a possibility that they 
may cause a virtual expansion of polygons due to these restrictions. The transformation process 
of the pivot move prohibits those moves in which self-intersections of polygonal segments occur. 
The angle 9 is selected from the range of —90° < 9 < 90°, excluding rotations with 9 > 90°. This 
tends to bring polygons to less folded conformations, and hence the polygon size is increased. 
Furthermore, the rotations are not completely randomly distributed over the range of 9, since 
the rotations are limited to those which make no change of the knot topology. 

In this paper, we modify the algorithm in two respects: We allow the selected chain of seg- 
ments to rotate by an angle 9 between and 360 degrees. We do not check the self-intersections 
during the process of rotation of the chain. In addition, we neglect the possibility of self- 
intersections occurring in the configuration after a pivot move is completed, since such configu- 
rations are negligible in the space C. 

With this algorithm, the topology of a polygon may change by a pivot move. We estimate 
the knot type of polygons by calculating several topological invariants |2.3| . the value of the 
Alexander polynomial Ax{t) at t = —1 [23], and the Vassiliev invariants of the second and 
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third order, v^K) and vs(K), respectively |25| . The chance of miss identification is negligible 
for the simple knots we are concerned with. 

We choose a certain conformation of a polygon as an initial state of the Markov process. The 
mean deflection angle between adjacent segments is equal to 90° in the thermal equilibrium. We 
prepare seeds on the cubic lattice as an equilibrium conformations of a polygon, exploiting the 
method used for the SAP on the cubic lattice [2B]. Starting from a seed, we generate a sequence 
of polygons by applying the pivot moves repeatedly. After discarding the initial 2000 transient 
conformations, we take samples of polygons at every 200 pivot moves. 

The simulation has been performed for polygons with the length N between 100 and 2200 
and for each of three knot types 0, 3i and 3ijJ4i. We have collected 10 5 polygons for each given 
length N and for each given knot type and we have evaluated the radii of gyration (J1J) and ©. 
We put a lower bound N m i n in fitting the simulation results using the scaling formula 1(H)). Thus 
we make the fit in the range N m i n < N < N max {= 2200), varying the value of N m i n from 100 
to 600. 

4 Results of The Simulation 

We now discuss the results of our simulation for the three knot types. The best fit curve for 
each knot type, together with the data on the ratio R 2 K /R 2 , are shown in Fig^ We have found 
that the scaling relation (jEJ) fits the simulation data very well for all three cases in the range of 
N from N m i n > 100 to N max = 2200. The estimated values of the parameters in @ are shown 
in Tabled The x 2 values remain small, x 2 — 1 _ 2 per datum, even if N m i n is lowered to 100. 
In the case of 3i, however, the x 2 values becomes small for N > 600. 

The ratios R 2 K /R 2 increase as functions of N, and they become larger than one for large N. 
The same behavior has been observed in other models, such as the Gaussian random polygon 
model |llj and the cylindrical SAP model |12j . The above observation tells that the average 
size of polygons grows due to the topological constraint. The exponent vr obtained in the fit 
are insensitive to the change of N m i n , showing the validity of the scaling relation R 2 K ~ N 2uk . 
The resulting values 2uk — 1-11 - 1.16 are consistent with vsaw '■> supporting the assertion in 
the papers [T§ll9*]. 
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We have found another optimal fit to our simulation data, as shown in Fig|21 The curves 
represent the second fit in which only data for N > 600 are used in Fig[2l The estimated values 
of the parameters are given in Table |2J The x 2 values for these fits grow rapidly as N m i n is 
lowered less than about 400. It implies that the formula © does not give a good approximation 
in the range N < 400. 

1.8 
1.6 
1.4 

~ 1.2 
1 

0.8 

100 1000 

N 

Figure 1: Simulation data on the ratios R K /R 2 for the knot types , 3i, 3ift4i and the best fit 
curves (solution 1) using ©. The scales are double-logarithmic. 

The second solution implies that the asymptotic behavior of the average size of polygons is 
scarcely affected by the topological constraint. It gives 2vk — 1.01 - 1.07, which is not much 
different from 2v = 1. Another noticeable difference is in the value of Ak/A. Ak/A < 1 for the 
solution 1, whereas Ak/A ~ 1 for the solution 2. 

If we assume Ak/A = 1, then the effect of topological constraint could only appear in the 
coefficients Bk- It is noted that the sign of ABk is negative for each knot type with large 
statistical errors. This observation supports the result derived by the perturbation argument in 
the limit of large characteristic length j2jj . 

To summarize, we have confirmed that the simulation data on the average size R K of ring 
polymers can be fitted by the scaling relation @. We have found two solutions. The best fit 
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Table 1: The value of the parameters of the scaling relation © of the ratio B? K / E? obtained 
from the best fit. 

trivial / all 





A*/A 


Ai?0 


2Ai/ 


x 2 


100 


0.471±0.014 


2.197±0.159 


0.133±0.004 


0.86 


200 


0.405±0.027 


3.172±0.429 


0.151±0.008 


0.69 


400 


0.376±0.079 


3.779±1.619 


0.159±0.024 


0.70 


600 


0.242±0.154 


7.812±6.457 


0.207±0.068 


0.679 


3\/all 




A, J A 


A5 3l 


2Az. 3l 


x 2 


100 


0.494±0.030 


0.258±0.269 


0.130±0.008 


6.26 


200 


0.455±0.064 


0.722±0.782 


0.140±0.017 


6.56 


400 


0.624±0.221 


-1.281±1.973 


0.104±0.0420 


6.54 


500 


0.578±0.262 


-1.041±2.683 


0.114±0.053 


5.15 


600 


0.793±0.277 


-2.417±2.022 


0.074±0.041 


1.30 


3itJ4i/aZZ 




A^aJA 


A5 3ltt4l 


2Az/ 3lB4l 


x 2 


100 


0.460±0.048 


-1.229±0.443 


0.140±0.013 


1.68 


200 


0.389±0.064 


-0.369±0.831 


0.160±0.020 


1.52 


400 


0.378±0.159 


-0.110±2.559 


0.163±0.049 


1.64 


600 


0.293±0.335 


1.670±8.634 


0.192±0.128 


1.46 
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Table 2: The value of the parameters of the scaling relation @ of the ratio B? K /B? obtained 
from the best fit. 



trivial / all 




A $ /A 


AB % 


2Ai/ 


x 2 


200 


1.720±0.464 


-3.731±0.915 


-0.029±0.035 


49.42 


400 


1.736±0.585 


-4.820±1.37 


-0.023±0.042 


7.24 


500 


1.359±0.270 


-5.345±0.844 


0.012±0.024 


1.91 


600 


1.203±0.514 


-3.884±2.214 


0.026±0.050 


1.66 


3i/ all 




AsJA 


AB 3l 


2Au 3l 


x 2 


200 


1.273±0.340 


-3.867±0.941 


0.012±0.034 


61.19 


400 


1.117±0.370 


-4.181±1.461 


0.034±0.040 


9.17 


500 


1.286±0.261 


-5.107±0.882 


0.019±0.025 


1.92 


600 


1.095±0.360 


-4.180±1.684 


0.037±0.039 


1.50 


3iJ)4i/a// 




As^JA 


A J B 3ltt4l 


2Az/ 3ltt4l 


x 2 


200 


0.926±0.306 


-4.056±1.137 


0.050±0.043 


12.22 


400 


0.997±0.432 


-4.992±1.787 


0.047±0.053 


3.80 


600 


0.877±0.827 


-5.021±4.623 


0.066±0.111 


2.67 
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1000 



Figure 2: Simulation data on the ratio B? K /R 2 for the knot types , 3i, 3ij}4i and the curves of 
the second optimal fit (solution 2) using (jfifl. The curves represent the second fit in which only 
data for N > 600 are used. 

solution (solution 1) is good for 100 < N < 2200 with quite small \ 2 (x 2 — 1) while the second 
solution (solution 2) becomes good for N m i n > 600. The two solutions lead to two different 
interpretations regarding the effect of fixed knot topology on the average size of ring polymers 
as their length N increases. 
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